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Abstract 

We establish the unimodality and the asymptotic strong unimodality of the ordinary 
multinomials and give their smallest mode leading to the expression of the maximal 
probability of convolution powers of the discrete uniform distribution. We conclude 
giving the generating functions of the sequence of generalized ordinary multinomials 
and for an extension of the sequence of maximal probabilities for convolution power 
of discrete uniform distribution. 

Since the eighteenth century, the expression of the convolution power of the discrete uni- 
form distribution has been very well known (e.g. de Moivre in 1711, see jTTJ 3rd ed., 1756] 
or [TJ1 1731]). This probability distribution arises in many practical situations including, 
in particular, games with equal chance, random affectation of tasks for many servers, and 
random walks. It is well known, see Dharmadhikari & Joak-Dev 7', 1988, p. 108-109.], that 
the convolution of two discrete unimodal distributions may be non unimodal. However, if 
these distributions are symmetric, we obtain a symmetric unimodal distribution. It is a 
discrete analog of Wintner's Theorem [191 1938]. Knowing that the convolution power of 
the discrete uniform distribution is symmetric unimodal, the determination of the maximal 
probability (mode) of such a distribution and its argument remains a question for consider- 
ation. As a recent work on the problem, one can see the article by Mattner & Roos [9, 2007] 
where they establish the upper bound for the maximal probability c 9i l < y/d/^q (<Z + 2) L 
(cq } L being the maximal probability of the L-th convolution power of the discrete uniform 
distribution on {0, 1, q}). Before them, there were several works aiming at finding such a 
bound. For example, Siegmund-Schultze & von Weizsacker [T5, 2007] proved the existence 
of a constant A such that c qi L < A/ (q + 1) \[L and gave an application of these upper 
bounds in the construction of a polygonal recurrence of a two-dimensional random walk. 
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Alternatively, our aim is to give an explicit expression of the mode of the L-th convolution 
powers of the discrete uniform distribution (section two) by means of the unimodality of 
the ordinary multinomials for which we study also the strong unimodality (section one), 
we end the paper (section three) by giving the generating functions for the two sequences 
of generalized ordinary multinomials: {(^) g }n and {(™ n z )J n , zgC, and thus of {c q ^n/q}n: 
when q is even. 



1 Unimodality of ordinary multinomials 

The ordinary multinomials are a natural extension of binomial coefficients (see [U 2007] 
for a recent overview on ordinary multinomials). Letting q, L € N, for an integer k = 
0, 1, . . . , qL, the ordinary multinomial is the coefficient of the fc-th term of the following 
multinomial expansion 

k>o 1 

with ( k ) 1 = (k) (being the usual binomial coefficient) and (^) — for k > qL. Using the 
classical binomial coefficient, one has 

0l+h-\ hj q =a 



(i + x + x 2 + --- + x«y = u ) xK (1) 



Readily established properties are the symmetry relation 

k) a = (qL-k 



(3) 



and the recurrence relation 



3 - £ (* -1) • w 

As an illustration of the latter recurrence relation, we give the triangles of pentanomial 
and hexanomial coefficients which are just an extension, well known in the combinatorial 
literature, of the standard Pascal triangle. 

Table 1: Triangle of pentanomial coefficients: (^) 4 

L\k 1 2 3 4 5 6 7 8 9 10 11 12 13 

1 

1 11111 

2123454 3 2 1 

3 1 3 6 10 15 18 19 18+ 15+ 10+ 6+ 3+ 1 

4 1 4 10 20 35 52 68 80 85 80 68 =52 35 20 ••• 

5 1 5 15 35 70 121 185 255 320 365 381 365 320 255 ••• 
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Table 2: Triangle of hexanomial coefficients: ( fc ) g 
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Let us investigate the unimodality of the sequence {(,) }^L - A finite sequence of real 
numbers {a/cjJJLo ( m — -0 ^ s cai l e d unimodal if there exists an integer I G {0, ...,m} 
such that the subsequence {ajt}j_o increases, while {afeltL; decreases. If a Q < a\ < ■ ■ ■ < 
a; _i < ai a = ■ ■ ■ = ai t > a^+i > • • • > a m then the integers Zo , . . . ,l\ are the modes of 
{ a fe}/JLo- I n * ne case wnere = h w< 3 talk about a peafc, otherwise the set of the values 
of the mode is called plateau. For positive non increasing and non decreasing sequences, 
unimodality is implied by log-concavity. A sequence {a,k}™ =0 is said to be logarithmically 
concave /log-concave for shortj or strongly unimodal if af > a;_ia;+i, 1 < I < m— 1. Also 
if the sequence is strictly log-concave (SLC for short), i.e. if the previous inequalities are 
strict, then the sequences have at most two consecutive modes (a peak or a plateau). For 
these notions, one can see Belbachir & Bencherif [3l 2007], Belbachir & al [2, 2007], Bertin 
& Theodorescu [I, 1984], Brenti [SJ 1994], Comtet 1970], Dharmadhikari & Joak-Dev 
[71 1988], Keilson & Gerber 1971], Medgyessy QH 1972], Sagan [14, 2007], Stanley [HI 
1986] and 16, 1989] and Tanny & Zuker PH 1974]. In the following, [o| denotes the 
greatest integer in a. 



The first main result of this article is the following. 

Theorem 1 Let q > 1 and L > be integers. Then the sequence {(j^) q }k=o * s un i' m 'Odal 
and its smallest mode is given by 

k L := argmax Q = [{qL + 1) /2J , 
Furthermore, we have the following recurrence relation 




where 

{{-q/2, . . .,q/2} ifq is even, 

{- (q + 1) /2, . . . , (q - 1) /2} if q and L are odd, 

{- (q - 1) /2, . . . , (q + 1) /2} otherwise. 

Proof. It suffices, for each one of the two cases: q odd and q even, to proceed by induction 
over L using the recurrence relation ([4]). □ 

Remark 2 For odd qL we have a plateau of two modes: qL/2 and qL/2 + 1. Otherwise 
we have a peak (qL + 1) /2. 
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2 Determining the maximal probability for convolution 
powers of discrete uniform distribution 

We are now able to achieve our purpose: the expression of the maximal probability of the L- 
th convolution powers of the discrete uniform distribution. Let U q be the random variable 
of the discrete uniform distribution on {0,1,..., q} and let U* L be its L-th convolution 
powers: 

1 

<?Tl 



U q := (S + Si + ■ ■ ■ + S q ) (S a is the Dirac measure). 



In [TJ 2007], Belbachir and al. established a link between the ordinary multinomials and the 
density probability of convolution powers of discrete uniform distribution. With respect to 
the counting measure, such a density is given by 



P(U* L =k)= KkJq fe = 0,l,...,gi. 
+ 



Remark 3 From Odlyzko and Richmond \1<A 1985] we know that for L sufficiently large, 
the sequence of probabilities {P (U* L = is strongly unimodal, from which we easily 

deduce that the sequence {(^) }k is also asymptotically strongly unimodal. 



Conjecture 4 For each positive integer q, the sequence {( k ) q }k is SLC. 

From Theorem^ as second main result, we give the values of c 9i l := max^ (^) / {q + 1) L ■ 

Theorem 5 The maximal probability of the L th convolution power of the discrete uniform 
distribution over {0, 1, . . . , q} is 

1 ( L 

Cq,L 



(q+l) L \l(qL + l)/2\ 



3 Some generating functions 

As a third main result, we give the generating functions for the sequence of generalized 
ordinary multinomials, the sequences {(^) q }n and {{"^) q }n, z E C, and the extended 
sequence of maximal probabilities for convolution power of discrete uniform distribution: 

Definition 6 For z G C, we define the generalized ordinary multinomials, as follows 

(A \p z(z-l)---Q-fci + l) f5 s 

W," , 1+ ^ +fe ^( fc 1 - fc 2 ) ! (fc 2 -fc3)!---(Vi-^)!fc 9 !' 

This definition is motivated by the relation (0). 
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Lemma 7 We have the following inequality 

z (z - 1) • • ■ (z - ki + 1) z(z- !)• • • (z- ki + 1) 



Ez yz — i) ■ ■ ■ yz — ki -i- L) _ >r-^ 



fcH hfc, = fc v ; V 9 ± 9/ 9 fti+2ha+— 9+9 

hi+h,2-{ hfc 9 =fei 

Theorem 8 Lei z£C, t/ie generating function for generalized ordinary multinomials is 
given by 

Wf) r = (i + * + * 2 + - + *T- 

n>0 ^ ' 9 

Proof. Using the Lemma, we have £„ >0 CO/ 1 = E (*) osrfes? *" 

h 1 +2h 2 -\ \-qh q =n 

h\+hz-\ Vh q =m 

On the other hand 

(i +t+t a+. ..+*,)* =E m >„e l )(t+^ 2 +---+^) m 

— l^ m >0\m) Ji!/2!-/»-il/,,I 

ii+2/ 2 H h9i«=™ 

ii+i2H \-lq=m 

We conclude by summation over n > is equivalent to summation over m > 0. □ 
Remark 9 Problem 19 of Comtet Jtjjjj, Vol.1, p. 172, states that 

J2 x n C t « (l + t + t 2 )' 1 = (1 - 2x - zx 2 y h , 
n>0 

using the fact that the coefficient oft n in the development o/fl + 1 + £ 2 )" : C f n (l + 1 + t 2 )™ = 
(™) 2 is maxfe (^) 2 , we obtain the following combinatorial identity 



-1/2 



-1/2 



3, 

T/iis lasi identity can be shown as the generating function of the sequence {c2, n }n- 
Theorem 10 Let z £ C, t/ie generating function of the sequence {( nz ) q }n is given by 

r nz\ ( u + 2u 2 + • • • + qu 9 



*rt\n ) . V l + ?i + tt 2 H \-ui 

n>0 9 

where u is a solution of the equation t = u (l + u + u 2 + ■ ■ • + u 9 ) z . 

Proof. Use Hermite's Theorem [6] for the function t >— ► t (1 + 1 + ■ • • + t q )~ z . □ 
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Theorem 11 For q even, the generating function of the sequence {c q 2 n /q}n is given by 

r (A.-TTtnn (i 2 ^ + ^ 2 + --- + g^ V 1 _g i + Et 2 i("" fc + " fc ) 

'k ql + u + u2 + --- + u*) 2 ' 

where u is a solution of the equation 



t = u 



9 + 1 \ 2/q ( g + i \ 2/? 



1 + u + u 2 H hu' 



Proof. Use the above Theorem for z = 2/q, and the change of variable t — > (q + 1) t. □ 
Remark 12 T/ie sequence {c q 2 n /q}n contains strictly the subsequence {c q i,}L- 
Corollary 13 For q = A, the generating function of { c^^ln is given for t £ ] — v5, 1 [ 

G 4 (t) := E t% n/a = (l - ^ 2 - ^ 4 - ^-t (5t 2 + 20) f ) V2 . 

n>0 ^ ' 

Corollary 14 PFe /iawe t/ie following identities 

(-5)- n f n/2 ) =2 and ^=2/^5. 



n>0 x ' 4 n>0 



Remark 15 TTie generating function of the sequence {c^ n } n is given for t G ]— 1, 1[ by 

n>0 
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